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Abstract 



A calculation of the inelastic scattering rate of Xe atoms on Cu(lll) is pre- 
sented. We focus in the regimes of low and intermediate velocities, where the 
energy loss is mainly associated to the excitation electron-hole pairs in the 
substrate. We consider trajectories parallel to the surface and restrict our- 
selves to the Van der Waals contribution. The decay rate is calculated within 
a self-energy formulation. The effect of the response function of the substrate 
is studied by comparing the results obtained with two different approaches: 
the Specular Reflection Model and the Random Phase Approximation. In 
the latter, the surface is described by a finite slab and the wave functions 
are obtained from a one-dimensional model potential that describes the main 
features of the surface electronic structure while correctly retains the image- 
like asymptotic behaviour. We have also studied the influence of the surface 
state on the calculation, finding that it represents around 50% of the total 
probability of electron-hole pairs excitation. 
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An incident particle scattering with a solid metal surface losses energy mainly through 
two channels: excitation of phonons and electronic excitations (electron-hole pairs, plasmons, 
etc) . At low energies, creation of phonons and electron- hole pairs becomes the main available 
channel. Understanding the different mechanisms for energy loss is a difficult task that has 
attracted both experimental and theoretical work. In an early experiment, Amirav and 
Cardillo JIJ reported a method for measuring directly the electron-hole (e-h) pair excitation 
in the scattering of Xe atoms on a semiconducting surface. More recently, it has been 
possible to measure the e-h pairs created by adsorption of hydrogen and deuterium 0-[4j] on 
transition metal surfaces. 

In the theoretical side, e-h pair creation probabilities due to moving atoms have been 
calculated with the use of several models. The most simple calculations were carried out 
within the Specular Reflection Model (SRM) ||, in which the surface potential is replaced 
by an infinite barrier ||. This implies that the electronic density vanishes at the surface, so 
the surface polarizability is not well described and energy losses are underestimated. More 
realistic calculations were performed by using the jellium model |7HTT| but surface states, 
which represent an important channel of electron-hole pair creation, are not included in 
this description. Very few calculations of the energy loss to date incorporate most of the 
complexity of the surface. One example, however, is the very recent ab initio calculation of 



the e-h pair creation by hydrogen scattering on Cu(lll) performed by Trail et al. [12]. There 



is also a large amount of theoretical calculations of energy loss of charged particles, ranging 



from very simplified models for the response function |1J-|T8| to more accurate self-consistent 
jellium calculations • 

In this paper, we calculate the contribution of the long-range van der Waals term to the 
probability of excitation of e-h pairs by an incident atom interacting with a metal surface. 



This problem was also studied by Annett and Echenique pi ], in an attempt to explain 
the experimental results of Amirav and Cardillo. They used a self-energy formalism, with 
the response function described within the SRM, and obtained a probability four orders 
of magnitude less than the experimental one. Several reasons might explain this large 
difference. On the one hand, the use of the simple SRM model for the response function. 
On the other hand, several other contributions to the scattering process were not taken into 
account, like those related to the atom-surface overlap (shorter-range effects), interaction 
between the static dipole generated in the atom by the surface proximity and the substrate, 
and non-adiabatic effects. We also focus on the van der Waals term, leaving the exploration 
of the other terms to a future publication, and try to study the influence of the model 
used for the response function, the semiclassical SRM or the more realistic Random Phase 
Approximation (RPA). RPA provides a more correct description of the electron density at 
the surface region and permits to investigate the importance of the surface state band in 
the electron-hole pair excitation. 

Unless otherwise is stated, atomic units are used throughout, i.e., m=e 2 =h=l 
We consider neutral atoms moving parallel to a metal surface with a constant, non rel- 
ativistic velocity v, at a distance Zq from the substrate, which occupies the z < half- 
space. We will start from the Born-Oppenheimer wave function for the atom, \a >= 
(j) a (R,)il} a (ri, T2, R), where the rj are the electronic coordinates and R the nuclear ones. 
Ground state wave functions of the solid will be denoted by \(3 >, and \f3' > will be the 
excited states. Neglecting coupling between the atom and the substrate, the many-body 
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wave function of the system can be written as the product of the respective wave functions 
of the atom and the solid. 

The contribution to the energy loss due to the long-range van der Waals coupling between 
the surface and the atom can be obtained by defining a self-energy, whose real part represents 
the van der Waals energy, and its imaginary part is related with the inelastic scattering rate, 
T, so that T = — 2Jm(S). The transition rate of the incident atom by excitation of electrons 
from the Fermi sea to unoccupied levels can be calculated by introducing the Coulomb 
operator, V, as the coupling between atom and substrate in Fermi's golden rule. Therefore, 
the self-energy will be given by: 

\<a(3\V\a'(3'>\ 2 
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Where the e are the corresponding energies of the states of the incoming atoms and the 
substrate. We can simplify eq.(|l]) by introducing the response function of the metal and the 
two-dimensional Fourier transform of the Coulomb potential, obtaining: 
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Where u is the transferred energy, and p a (r) is the charge-density operator for the atom. 
Eq.(^) is valid for arbitrary trajectories. Now we focus on trajectories parallel to the surface, 
replace the density operator by p a (r) = Z5(r — R) — J2i^( r — r i) ; where Z is the nuclear 
charge, and take nuclear wavefunctions as, 

0a(R) = j^y^ {tMvn) ^ - *) (3) 

Assuming the atom is far from the surface, we can make a multipole expansion in which 
we take the dipole term as the leading contribution, because for neutral atoms there is no 
monopole contribution. This term introduces the dipole oscillator strengths, corresponding 
to transitions from electronic state to n, f n0 = 2u n0 \ < n\r\0 > | 2 /3 where u n Q are the 
frequencies of the dipolar transitions. These transformations in eq.(0) lead to the following 
expression of the scattering rate: 
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(1 + sig(z - z')sig(z - z"'))5{u + u n0 - q||-V|| + q 2 /2M) 



Where sig(z-z') stands for the sign function. Once the energy transfer is fixed by 
Dirac's delta, maximum and minimum momentum transfers are given by = M{vcosO ± 



v 2 cos 2 8 — 2uj n o/M). 6 is the angle between the parallel momentum and the initial velocity. 
M is the mass of the atom, which for Xe atoms is of the order of 10 4 a.u. Energies above 
the 9 eV threshold, corresponding to the first available dipole transition in the Xe atom, 
are needed so that e-h pair excitation is feasible. Therefore, it is the gap in the atomic 
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levels which determines the need of a minimum momentum transfer different from zero. q + 
corresponds to the transfer of the whole initial energy. 

We have used two different models for the response function. In the SRM the surface 
potential is replaced by an infinite barrier and it is assumed that the electrons are specularly 
reflected at the surface, with no interference terms. Therefore, the surface response function 
is strictly zero beyond the jellium edge. In the more realistic RPA model [f23H 25|, electrons 
respond as a non-interacting gas to the total potential, and the response function is obtained 
solving a self-consistent equation on the induced potential. We model the metal by a slab, 
such that one-particle wave functions in the directions parallel to the surface are plane 
waves, while in the perpendicular direction they are the solutions of a one-dimensional 
model potential [p6|| . This potential reproduces the width and position of the band gap, 
as well as the binding energies of the surface and image states. Close to the Fermi level, 
the resulting wave functions are comparable to those obtained with an ab initio calculation 

. The potential has also been proved to accurately predict the broadening of surface and 
image states on metal surfaces p7|- p0|1 . 

For the RPA calculations we have used a supercell geometry where the Cu(lll) surface 
is represented by a 30-layers slab, and the vacuum region is equivalent to 20 atomic layers. 
In the SRM model M, the surface response is directly obtained from Mermin's dielectric 
function [p2]| . We have used the values 2.67 and 0.013 a.u.for r s and the damping parameter 
7 respectively. Our results exhibit only a slight dependence on 7. This is a clear indication 
that, for the moderate projectile velocities considered here, the main contribution to the 
inelastic scattering rate (r) comes from the creation of e-h pairs. Therefore, from a practical 
point of view, it is justified to consider V as equivalent to the probability of e-h pair creation, 
and we will make no distinction in what follows. The energies and oscillator strengths of 



the electronic transition of the Xe atom were taken from Ref. 



Fig.p] depicts V for trajectories parallel to the metal surface with different incident veloc- 
ities as a function of the impact parameter zq (measured from the jellium edge) calculated 
using a) the SRM and, b) our model RPA surface response. The e-h pair creation probability 
obtained within the SRM approximation is roughly three times smaller, and decays much 
faster with zq than the one calculated using the RPA model response. These results indicate 
that the polarizability of the copper surface is seriously underestimated in the simplified 
SRM model due to the presence of a fictitious infinite barrier at the surface. In fact, the 
induced electronic charge becomes strictly zero beyond the jellium edge. This is in contrast 
with the model RPA response, which peaks at the surface and penetrates some atomic units 
into the vacuum (see Fig.[3]), allowing for higher scattering rates for larger impact parame- 
ters. Our data also point out that, in order to obtain a reliable description of the inelastic 
scattering at low and intermediate velocities, it is necessary to explicitly include a relatively 
accurate description of the electronic structure in the vicinity of the Fermi level. This is 
specially true for the surface band, as will become clearer below. 

In Fig.|2| we show the scattering rate for both SRM and RPA response functions as a 
function of the velocity of the incident atom, for two different impact parameters close to 
the surface. At very low v, no creation processes are possible because the initial energy must 
at least equal the energy of the first allowed transition in the incoming atom (~ 9 eV for Xe). 
When v increases, this threshold can be overcome, but there is a simultaneous requirement 
of small energy and relatively large momentum transfer to the target that severely reduces 
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the phase space of the allowed electronic excitations in the substrate. Increasing v the e-h 
pair creation probability increases until it reaches a maximum for v around 2 a.u.. The 
position of this maximum is only slightly dependent on the impact parameter. By further 
increasing v, T starts to decrease since, as the energy transfer gets larger, the low q\\ region 
available in phase space becomes smaller, and the surface response tends to decrease for 
large values of q\\. Finally, in the limit of very high v (not reached here) only the plasmon 
contribution becomes important. 

The surface state in Cu(lll) crosses the Fermi level, so it is a partially occupied band. 
This, in addition to the larger weight of this state in the surface region, makes reasonable 
the expectation that a large fraction of the electrons excited by incoming Xe atoms will 
be originated from inter and mainly intraband transitions in this band. We will study the 
effect of the surface band in the following. Fig.[3] shows the imaginary part of the RPA 
response function for our slab model Imx RPA (q\\,w, z, z' = z) for gn=0.3 a.u., uj—1 eV, 
as a function of z (coordinate normal to the surface), with and without the surface state 
contribution. The imaginary part of the RPA response function gives the number of available 
electronic transitions in the substrate. While the removal of the surface state does not affect 
the response in the bulk region, the number of excitations (e-h pairs in this energy range) 
generated at the surface decreases dramatically. 

Fig.[| shows T versus the impact parameter with and without the surface state, for an 
incident velocity of 2 a.u. Close to the surface, the result without the surface state is reduced 
to a 50% of what we get with all the states. This reduction, however, gets somewhat lower 
as the atom is farther away from the surface. 

In summary, we have reported calculations of the inelastic scattering rate of neutral 
atoms with moderate velocities by metal surfaces. We have restricted here to the long- 
range van der Waals term. Results for other contributions to the total inelastic rate will 
be presented elsewhere. We have found that the introduction of a more realistic model of 
the surface, which explicitly includes the electronic structure of the surface near the Fermi 
energy and where the response function is calculated within RPA approximation, increases 
the scattering rate given by the semiclassical SRM model by approximately a factor of 
three. These differences are enhanced at low and intermediate velocities. By removing 
the contribution of the surface state in the calculation of the RPA response function, we 
have seen that this band is the main source electronic excitations in the surface region, 
we conclude then that surface bands have to be correctly described in order to achieve a 
complete description of surface-atom scattering events. 
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FIGURES 



FIG. 1. Inelastic scattering rate for Xe atoms moving parallel to a Cu(lll) surface as a function 
of the impact parameter. The incident velocity is v=l a.u. (dot-dashed line), 2 a.u. (dashed) and 
v=3 a.u. (solid), a) results for the SRM surface response function, and b) for a model RPA 
response. 

FIG. 2. Inelastic scattering rate for Xe atoms moving parallel to a Cu(lll) surface as a function 
of the incident velocity. The impact parameter zq is 1 a.u. (solid lines) and 2 a.u. (dashed), a) 
results for the SRM surface response function, and b) for a model RPA response. 

FIG. 3. Imaginary part of our model RPA response function for q\\=0.3 a.u. and u=l eV, as a 
function of z (coordinate normal to the surface), for z = z' with (solid line), and without (dashed 
line) the surface state contribution. 

FIG. 4. Inelastic scattering rate for Xe atoms moving parallel to a Cu(lll) surface as a function 
of the impact parameter. The model RPA surface response function was used with (solid lines) 
and, without (dashed) the contribution of the surface state. 
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